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states at arbitrary filhng fraction possess a degeneracy with respect to center 
of mass translations. An analysis is carried out for multi-component elec- 
0^ . tron systems and extra degeneracies are shown to exist. Their application to 
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numerical studies is discussed. 



I. INTRODUCTION 

Following Lauglilin's original work [jT|, Halperin 0| proposed an extension to the orig- 
inal Lauglilin-Jastrow wavefunction incorporating the possibility that the electrons have 
some further degree of freedom in addition to their two-dimensional coordinate, which in 
Halperin's original observation was taken to be spin. Recent experiments [^ on electron 
systems where the extra degree of freedom can be seen to indicate the state of the electron 
in the third direction seem to suggest the existence of new universality classes of states, 
as first proposed by Haldane and Rezayi [Q and analyzed by Yoshioka, MacDonald, and 
Girvin [^. The possibilities raised by this extra degree of freedom have been analyzed in 
various contexts [§-|l9|. 

In the standard treatment of the problem of electrons in a magnetic field subject to 
periodic boundary conditions (PBCs), Haldane has shown that every eigenstate of a 
translationally invariant Hamiltonian has a degeneracy with respect to center of mass trans- 
lations. For a state at filling fraction p/q the degeneracy is simply q fold. Haldane and 
Rezayi further constructed the explicit generalization of Laughlin's wavefunction to the 
PBCs, again showing the q fold degeneracy inherent to states in this geometry. In light of 
the interesting possibilities raised by the double layer systems, it is natural to generalize 
the wavefunction construction to multi-component systems subject to periodic boundary 
conditions. 

The remainder of the paper is organized as follows. In the next section we review the 
construction of the single layer generalization of Laughlin's wavefunction subject to peri- 
odic boundary conditions. We then extend this construction to multi-component systems, 
showing interesting new degeneracies distinctive to these systems. Finally, we discuss the 
quantum numbers of these states with regard to numerical studies and their use in distin- 
guishing between possible ground states. 



II. LAUGHLIN WAVEFUNCTION IN PBC 

We shall construct the Laughlin state subject to periodic boundary conditions (PBC's), 
closely following [Q. In a system subject to PBCs in a magnetic field, the translation 
operator takes the form 

t(a)=exp(a-(V-^A)-z^) (1) 

and obeys the non-commutative algebra 

t(a)t(a') = t(a + a')e'^ (2) 



where / = J{fi/eB) is the magnetic length. We wish to impose generalized boundary 
conditions by requiring that all physical quantities be invariant under translation of any 
particle by the set of translations L^n = mLi + nL2 where 

|Li X L2I =27rA^0/2 (3) 

and N(i) is the number of flux quanta. We impose the general boundary conditions on the 
wavefunction for any particle i 

t,(L^J^ = (r/^O'^^e**"''''""^ (4) 

where rfmn = (— 1)(™'+"+™'^) and we will choose $0 = as our boundary condition. The 
physical region under consideration can be seen to be deflned by four points z = |Li(±l±r) 
where r = j^e'^^. We shall use the symmetric gauge A = (—By/2, Bx/2). The single particle 
wavefunctions in the lowest Landau level are then given by 

ip{x,y) = e~^^^^ f{z) z = X + iy (5) 

where f{z) is an analytic function with N^ zeros in the principal region. Applying the 
boundary conditions to (|) we flnd 

f{z + L^n) _ (^^ j^^exp f ^^"^"^" ,"t ^^"^""' 'l (6) 



f{z) """'' "V 4/ 



where we write a two-dimensional vector a as a^ + lay. The basic building block that we 
shall use is the quasi periodic function w{z) which obeys 

w{z) -^--^^P^ m^p >■ ^^^ 

An explicit representation of the function ^(2;) is given by 

w{z) = exp(^^^jei{Kz\T) (8) 

where Lmn = K~^7r(m+r?,r) and Qi{k,z\t) is the odd elliptic theta function. We can therefore 
construct the single particle wavefunctions 

^W}(r) = ijlv{r- a„)exp(^^)) (9) 



^Q = l 



where 






and 



(P(t) = w(2;)exp 



/ z*z 



(11) 



The number of linearly independent solutions of @ and can be discerned in the following 
way. Inserting the above form into the boundary conditions yields the constraint 

where Lp^ is restricted to be a primitive translation. There are therefore (N^)"^ possible values 
for a. One way in which to resolve the degeneracy is to form a superposition of states where 
every zero is shifted by the same amount in some primitive direction. This is equivalent to 
constraining the wavefunction to be invariant under a translation of the electron coordinate 
in this direction. Given one set of zeros {aa} which satisfies the boundary condition (|12|), 
another set can be generated by shifting the zeros 

tta^ aa + ^ (13) 



where L' is also a primitive translation. Let us then form a linear superposition of states 



pr 

with shifted zeros 



N4, 
vl/(r) = ^M/(r)|,„,^} (14) 

7=1 

where 

where Ki is a primitive translation. We can operate on the above wavefunction by perform- 
ing an overall shift of the zeros 

V7 = "«'7 + ^fc^ (16) 

where K2 is a primitive translation obeying 

|Ki XK2I = 27rA/'/. (17) 

We can therefore see that n^ = 1,...,A'^ using (P^). The number of linearly independent 
solutions is given by A^^. It is important to note that this is one program which generates a 
space of linearly independent solutions but not the only one. 

We now consider the many particle Laughlin wavefunction at Landau level filling v = 
p/q = 1/m on the plane 



i<j i 



4/2 



(18) 



In the following, we shall denote N^ = pN and A^,^ = qN. When the system is subjected to 
PBCs, this wavefunction generalizes to 

^(M) = (n¥'(ri-rj)") 

X n(^(R-^a„)exp(|^)) (19) 

where 



R=Eri (20) 

i 

and (y5(r) is defined in (pJ]). If we apply tlie translation operator to any of the coordinates 
and apply the boundary conditions, we find 



2— Id "-a -L'pr 



(21) 



ni m 

where Lp,. is restricted to be a primitive translation. The center of mass portion of the 
wavefunction can be seen to be formally equivalent to the single particle wavefunction for a 
particle of charge eNe which sees m flux quanta. By analogy with the single particle case, 
we therefore conclude that the Laughlin state on the torus has an m fold degeneracy related 
to the action of the center of mass translation operator. Very generally, this degeneracy is 
inherent to every eigenstate, as was shown by Haldane 0. 

III. EXTENSION TO MULTI-COMPONENT SYSTEMS 

We can write the planar extension of Laughlin's original wavefunction to a system with 
A^s species of electrons as first expressed by Halperin P| in the two-component case 

^^[{^.}] = X{{z^- z,r^-^'''^)X{eM-^,\z^?) (22) 

where ctj is the species quantum number and K is an N^ x A'"s symmetric matrix encoding the 
correlations between the electrons where we impose "DeiK > 0. If we consider the plasma 
analogy, we find that in order to have a uniform fluid stable against fluctuations, we must 
choose DetK > 0. We shall consider the case of DetK = separately. The matrix K also 
makes an appearance in effective theories of the fractional quantum Hall effect, as the long 
distance physics of the Hall fluid can be described by the Lagrangian 

^ = ;^(E^(^'^>''^'«M^X' + 2e'^'^'E^'^^A«A) 



Arc 



cr,cr 



+ Maxwell terms (23) 

where a'^ are A^^ Chern-Simons gauge flelds. This Lagrangian has been discussed extensively 
previously |jl3|. The two-component case is of special physical interest, where the two species 



of electron may represent the two possible physical spin states of an electron or, in double 
layer systems, which of two layers an electron resides in. We have suppressed the quasi-spin 
portion of the wavefunction, which in the case of fermions would correctly anti-symmetrize 
the overall wavefunction. 

We write the Halperin mult i- component wavefunction subject to PBCs as 

a ^ i 

" n(,-p(' — 1^^ — )j (24) 

where {S"^{a)} are integers. We also note the constraint that each electron must have N^ 
zeros in the wavefunction, implying 

N^ = Y.K{a,,a,) (25) 

3 

which correctly gives us the filling fraction 

y = Y.K-\c,,c,'). (26) 

We can apply the PBCs to any of the electrons to get the constraints on the center of mass 
portion of the wavefunction. Applying the translation operators to any particle and using 
the boundary conditions, we find 

m 

K{a, a') = ^ 5"(a)5"(a') m > N,. (27) 

a=l 

and 

^r,^"(a,)=L„,„(a,) (28) 

a 

where Lmni<^) is a primitive translation. Therefore, given that S"'{(j) is integer valued, this 
implies that the matrix K{a, a') must be positive definite, as is necessary for thermodynamic 
stability. We can also impose the constraint on the lowest common denominator 
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l.c.d{^"(ai), S''{a2), ..., ^"(cT^vJ} = 1 (29) 

for all values of a. 

We now will determine the number of linearly independent solutions to these equations, 
or equivalently, the number of independent ways of arranging the zeros of the center of mass 
wavefunction. For the purposes of forming sets of zeros such that the corresponding wave- 
functions are orthogonal, we may constrain the center of mass portion of the wavefunction to 
be invariant under independent translations of the individual centers of mass of each species, 
keeping in mind that we are not constraining the overall wavefunction to be invariant under 
these translations, only the center of mass wavefunction independent of the relative piece. 
In order to achieve this, we will form a superposition of states with shifted zeros as we did 
in the single particle case. We then wish to determine the set of allowed translations {7^} 
of the zeros such that 

r.-^r„ + 5:^"(a,)7, (30) 

j 
is allowed by the boundary conditions. Inserting this into the boundary conditions, we find 

7i = Yl K'^i^i^ (yk)Lmn{(yk)- (31) 

k 

By analogy with the single particle Hilbert space, the number of linearly independent solu- 
tions is simply given by 

Deg = DetK (32) 

As an example we consider the two component case where the K matrix is given by 

(mi n \ 
(33) 
n 1712 J 

and the degeneracy is mim2 — n^- This result was noted in the context of Chern-Simons 

effective field theory approaches to the Hall effect |T2|JT^ but not elaborated upon. The 



explicit result is crucial in numerical simulations, if one is to try to distinguish between 
possible ground states. The degeneracy provides another quantum number in addition to 
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filling fraction to determine the universality class of a given incompressible quantum fluid. 
It is important to make several comments. 

1. Every eigenstate of a Hamiltonian that depends only on the relative position of the 
particles has a degeneracy of q on the torus if the filling fraction is given by p/q. The 
degeneracy described here is only a feature of the state given by the Laughlin-Halperin 



form (^). The overall multi-component degeneracy is due to the independent translation of 
the zeros of the centers-of-mass of different species, subject to constraints from the correlated 
part of the wavefunction. This includes q translations of the overall center of mass of the 
system, as well as translations of the zeros which can be interpreted as relative translations 
of the centers-of-mass of the different species. 

2. We have assumed in the derivation that deiK > 0. If detK = and K{a, a') = m 
the system is equivalent to a single layer u = 1/m Laughlin state with an extra degree of 
freedom. In this case, the state has the usual center of mass degeneracy of m. In addition, 
the system has a residual U{1) degeneracy with respect to the number of particles of each 
species. While the total number of particles is conserved, the Hamiltonian is invariant under 
rearrangements of the number of particles of each species. This wavefunction has many other 



interesting features, as has been investigated in several papers |Ty,|T5|JT^. 

IV. QUANTUM NUMBERS OF MULTI-COMPONENT STATES 

In his analysis of two-dimensional electron systems in magnetic fields subject to PBCs, 
Haldane constructed a general symmetry formalism which clarified the center-of-mass 
degeneracy seen previously in finite size studies subject to PBCs. This construction has 
proven to be very useful in finite size studies on the torus, simplifying the calculations as well 
as providing a correct classification of states allowing a direct comparison with studies done 
in other geometries. We shall briefly review this analysis and proceed with its application 
to the multi-component systems. 

As before we shall consider a physical region defined by the vectors Li and L2 and require 
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physical quantities to be invariant under translation through L^n where L^n = ■mLi +nL2. 
We can consider the operator that generates relative translations between the particles while 
maintaining the center of mass, defined to be 

t.(pL_) = nt.(^)t,(^§=) (34) 

j e e 

where Nf. is the number of electrons in the system and the filling fraction is given by z/ = 
Ne/N^ = p/q and the operators tj(a) are the translation operators acting on particle i 
defined as before (p. We can rewrite the above operator, using the fact that the physical 
states that this operator acts upon are subject to the conditions (|) as 

ii{pLmn) = ti{pLmn)'[[tj{ TT^) 

= {Vmuf'^expitpi^o ■ L„„))T(-^) (35) 

where 

T(a) = nt.(a) (36) 

i 

is the center of mass operator. We can therefore classify the eigenstates of a translationally 
invariant Hamiltonian obeying [if, T(a)] = by the quantum number k defined to be 

Ti^m = (r^^„)^'^exp( /'^ + ^-|°^-^- )|v^). (37) 

As this operator commutes with the center of mass translation operator T(^^), each eigen- 
state of a translationally invariant Hamiltonian, labeled by k, will have q fold center of mass 
degeneracy. In the following we shall specialize to the boundary conditions with $o = 0. 

We can explicitly define eigenstates for the above operator, useful for numerical work, in 
the following way. We can denote a fundamental set of translations as 

ti^) = t, t(^J = t. (38) 

and using the relation 

2w 



tit2 = t2tie ^^ (39) 

10 



we can define a basis as 

ti\m) = e " ^<i> \m) t2\m) = \m + 1) (40) 



_ • 277 771 



where \m + A'^^) = \m). This is the smallest set of translations consistent with a given set of 
boundary conditions. If we denote the many particle occupation number state as 

\ni,n2,...,nNj (41) 

we can write the eigenstate of (|37|) in the following fashion 

1^) = Xll'^i + qk,n2 + qk,...)e^^''. (42) 



k 



If we act with the operator ( ^7\) on this state, we find 



T{-^)\\if) = (_i)'""P'?e-¥(J^-™+j'!''^)|^) (43) 



where 



Af0 

j. = Y.^ft. (44) 

t=i 

Therefore, to set up the Hilbert space, we group sets of states into their j^ value. We 
then construct the full Hilbert space by forming linear superpositions of these states, with 
their occupation numbers shifted, multiplied by a phase factor. These are our fundamental 
basis states. We note that j^ and jy are integers defined modulo N. If we then define two 
reciprocal lattice vectors as 

Gc,-Lp = 2tt5^p (45) 

one can see that the Brillouin zone pseudomomentum k is given by 

k=-(j.Gi+j,G2) (46) 

where 

Ux,3y) = {jx,Jy) - (y' y) P^ = odd 

= i3x,jy) pq = even (47) 
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The most symmetrical Bravais lattice that we can consider is the hexagonal PBCs, 
invariant under a 7r/3 rotation. We have defined the k = point in ( p?]) to be the unique 
point invariant under a vr/S rotation, where Li — > L2 and L2 —>■ L2 ~ Li. Therefore, the 
invariant point must be the solution to the equation 

r(^,.T(4)^r(|,^r(^,. (48) 

If we separate the last term and solve, we conclude 

T{^)\^{k = 0}) = {r^mnr\'^{\^ = 0}). (49) 

This defines the k = point. The significance of these comments for numerical work is in 
the fact that in the thermodynamic limit the k = point becomes rotationally invariant 
and therefore the signal for a QHE fiuid is a k = ground state with an energy gap to any 
excited states. 

We now turn our attention to the multi-component system. In this case, the same analysis 
carries through as above with the occupation basis being expanded to include the pseudospin 
quantum number. If one acts with the operator ( |5^ on the single layer wavefunction at 
u = 1/m, one finds the Laughlin state to be a k = eigenstate. In the following we shall 
denote DetK = qN'. If one acts with the operator (^) on the wavefunction \I'^[{2;j}] one 
can determine the quantum numbers 

j^=pNY,^p, (modiV) (50) 

i ^ 

and 

- A^- 
h=pNT.j;^P'. ("^odiV) (51) 

i ^ 

where pi and p'^ are specified integers and Ni is the number of electrons of species i. One 
can therefore see that in the multi-component systems, the Laughlin-Halperin like ground 
state will have an overall degeneracy given by Det-ft' with a center of mass degeneracy given 
by q as well as a A^' fold degeneracy related to the different translations of the centers of 
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mass of the different species. There are four points in k space that have parity invariance, 
{jx,jy} = {0, 0}, which corresponds to k = and three states that he on the zone boundaries 
{N/2,0}, {0, N/2} and {N/2,N/2}. The relative portion of the wavefunction is an eigen- 
state of parity so the k vector of the states representing a Laughhn-Halperin hke ground 
state can only lie on one of these four points. We can therefore write 

DetK = q{No + 3Nb) (52) 

where A^o is the number of k = states and Nb is the number of zone boundary multiplets. 
There are only two possibilities 

even integer Nq = Nb = N' /A (53) 



NiP 
for some species i, otherwise 



Nq = N', Nb = 0. (54) 



As an example, we examine the (3, 3, 1) state, represented by the matrix 

recently shown [^ to be a physically realizable state. In this case, as the filling fraction of 
the state is z/ = 1/2, we have a two fold center of mass degeneracy. From the above analysis, 
we also find that there should be degenerate states possessing j^,. = 0, N/2 and jy = 0, N/2. 
There are then four degenerate states corresponding to {jx,jy} = {0,0}, a k = state, 
{N/2,N/2}, {N/2,0}, and {0,iV/2} for an overall degeneracy of eight as expected. This 
result has been borne out by numerical studies |1^ . It should be noted that this degeneracy 



is a feature of the exact \E'''^[{2;j}] state, which can be generated numerically by choosing a 
truncated Hamiltonian possessing the correct pseudopotentials to make the \E'-'^[{zj}] state 
the unique ground state. 
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V. CONCLUSION 

In this paper it has been shown that the degeneracy of a multi-component quantum Hall 
state on the torus, denoted in the standard fashion by the matrix K is given by DeiK. The 
quantum numbers of these multi-component states in Haldane's symmetry analysis have 
also been determined. These predicted quantum numbers provide a powerful topological 
invariant with which one can distinguish possible ground states. One example of this is 
in the v = 5/2 system. There have been two wavefunctions suggested to explain the ob- 
served anomaly [^,^, both possessing the correct filling fraction. On the torus, however, 
in addition to their center of mass degeneracies, the Pfaffian state possesses a three fold 
degeneracy |2^ while the spin-singlet state has a five fold degeneracy |^. While this issue 



is not yet resolved, it is clear that the ground state degeneracy is a useful characteristic in 
distinguishing between these two states numerically. Lastly, there is a remarkable transition 
in the z/ = 2/3 double layer system between two different ground states that occur at the 
same filling fraction but with different degeneracies on the torus, as seen in numerical stud- 



ies 



20| . It provides a fascinating example of the importance of degeneracy considerations in 
numerical studies on the torus. 
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